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,—1- \ Abstract: We constrain the rotation angle a of the Hnear polarization of CMB photons 

•^ \ using the large angular scale (up to ~ 4°) signal in WMAP 7 year data. At these scales, 

the CMB polarization pattern probes mostly the reionization era. A genuine rotation 
f— ^ ■ may be interpreted as cosmological birefringence, which is a well known tracer of new 

^"^ ■ physics, through the breakdown of fundamental symmetries. Our analysis provides a = 

Tl ! -1.6° ± 1.7° (3.4°) at 68% (95%) C.L. for the multipoles range A£ = 2 - 47 not including 

.y—{. an instrumental systematic uncertainty that the WMAP team estimates at ±1.5°. This 

kS ! improves the bound obtained by WMAP team (Komatsu et al., 2010). Moreover we show, 

C^ \ for the first time at low multipoles, the angular power spectrum a£ in search of a possible 

scale dependence of the birefringence effect. Our findings are compatible with no detection 

at all angular scales probed here. We finally forecast the capabilities of Planck in tightening 

the present constraints. 
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1. Introduction 

The Cosmic Microwave Background (CMB) can be used to constrain fundamental symme- 
try breaking. Under the Standard Model of particle physics, the CMB temperature to B 
mode polarization correlation (TB) and the E mode to B mode polarization correlation 
(EB) are expected to vanish since parity is conserved [1, 2]. The latter prediction follows 
from the purely electromagnetic, charge blind nature of the physics of CMB anisotropics 
[3, 4] and becomes observable due to the opposite handedness of {T,E) vs B modes [5]. 
A non null measurement of either TB or EB would thus hint at new physics beyond the 
Standard Model (or suggest a careful reanalysis of the dataset in search for un-removed 
systematic effects). 

A number of models predict non-standard parity (P) and CP violations (where C stands 
for charge conjugation), as well as CPT violations (T standing for time reversal) and the 
associated breakdown of Lorentz invariance [6]. Parity violations in the electromagnetic 
sector are associated with an in vacuo rotation of the linear polarization direction of a 
photon [3, 4], an effect known as cosmological birefringence. Several tests have been set 
forth, either in earthly and orbital laboratories [7] or through cosmological observations 
[3, 4, 8, 9]. These violations should also impact the CMB, whose statistical properties are 
derived under the assumption of symmetry conservation. The CMB is a powerful probe 
for such analyses for two main reasons. First, it is generated in the early universe, when 



the physics at the stake was not obviously identical to present. Secondly, the long look- 
back time of CMB photons may render tiny violations to the electromagnetic Lagrangian 
observable, since such effects usually accumulate during propagation. 

CMB polarization arises at two distinct cosmological times: the recombination epoch 
{z ~ 1100) and the reionization era at z < 11 [10]. When the CMB field is expanded in 
spherical harmonics, the first signal mostly shows up at high multipoles, since polariza- 
tion is generated through a causal process and the Hubble horizon at last scattering only 
subtends a degree sized angle. The later reionization of the cosmic fluid at lower redshift 
impacts the low i instead. These two regimes need to be taken into account when probing 
for cosmological birefringence, since they can be ascribed to different epochs and, hence, 
physical conditions. 

Historically, the effect has been first constrained by measuring polarized light from 
high redshift radio galaxies and quasars [3, 4, 9, 11, 12, 13, 14], see [15] for a new analysis 
on ultra-violet polarization of distant radio galaxies. For recent, polarization oriented CMB 
observations [16, 17, 18, 19] it has become feasible to measure TB and EB correlations, 
other than TT, TE and EE correlations. While no detection has been claimed to date, 
polarization data have been used to derive constraints on the birefringence angle [17, 20, 
21, 22]. 

The goal of this paper is to constrain the cosmological birefringence angle at large 
angular scale considering the latest WMAP data release. Our method naturally provides 
also a scale dependence for such an angle. 

In the current analysis we do not address the issue of the systematic effects whose 
detailed analysis goes beyond the scope of the paper. Our new results are based on the 
WMAP data and we do rely on their original estimate of the uncertainty on the birefrin- 
gence angle coming from systematic effects [17]. For what concerns the Planck forecast 
we provide, we cannot give an estimate on the errors coming from Planck detectors since 
those are not yet publicly available. However, following [23], we will provide an estimate of 
the level of accuracy that polarimeter calibration has to satisfy in order to make the error 
coming from this systematic effect smaller than the statistical one. 

In the limit of constant birefringence angle, a, the angular power spectra of CMB 
anisotropics, assuming Cj^ = Cf^ = 0, are given by [1, 20, 29, 30] ^, 



(Cf^'°^^) = (Cr)cos(2a), (1.1) 

(Cj^'°^^) = (CD sin(2a) , (1.2) 

l^CfE,obs^ = (Cf ^) cos2(2a) + {Cf'') sin2(2a) , (1.3) 

l^CfB,obs^ = (Cf ^) cos2(2a) + (Cf ^) sin2(2a) , (1.4) 

(Cf ^'°^^) = \ {{Cn + (CD) sin(4a) . (1.5) 



^See [24, 25] as an example of computation that takes into account the time dependence of a in a specific 
model of pseudoscalar fields coupled to photons. See [26, 27, 28] as examples of non-isotropic birefringence 
effect. 



These equations follow from 

T.obs T , n,T /-, a\ 

«£m = «^m + «£m ' (1-6) 

o-fm^' = «fen cos(2q) + af^ sm(2a) + a";f , (1.7) 

«f^°'" = «ftn cos(2q) + af„ sm(2a) + a^^ , (1.8) 

where we have explicitly written the contribution due to residual instrumental noise present 
on the maps. As discussed above, the amount of rotation depends on the effective travel 
time of the photon and the effect is scale (and, hence, harmonic multipole i) dependent 
[31]. In the present paper, we focus on the low i contribution, with a twofold aim: from 
the one side, we provide a new constraint on a for low i using, for the first time (to our 
knowledge) a linear estimator in combination with an optimal angular power spectrum 
estimator (APS) for such a regime ^; from the other side, we give the first measurement of 
the spectrum of a at low i. Note that the technique used by WMAP team does not allow 
for possible scale dependence of a. 

Our constraint reads a = -1.6° ± 1.7° (3.4°) at 68% (95%) C.L. for Ai = 2 - 47. 
Considering A^ = 2 - 23 we obtain a = -3.0°t|fo at 68% C.L. and a = -3.0°t|fo at 
95% C.L.. This is the same multipole range considered by the WMAP team in [17] (the 
only other result available in the literature at these large angular scales) where with a pixel 
based likelihood analysis they obtain ^WMAPTyr ^ _3_go ^ 52° at 68% C.L.. 

The paper is organized as follows: in Section 2 we present the methodology adopted 
and describe the simulations performed and the tests proposed; in Section 3 we show the 
results obtained (i.e. a new and more stringent constraint of a and its spectrum at low 
i). In Section 5 we draw our main conclusions. In Appendix A further details about the 
behavior of the considered estimators are given. 



2. Methodology 

2.1 Estimators 

Assuming that the amplitude of the primordial B-mode is zero, we write down the following 
linear estimators [22] 

Dg = C^ cos(2a) — Cf sm(2Q) , (2.1) 

Dg^ = Cg ' - - I C^ + C^ ' I sm(4a;) . [2.2) 

The linearity of the expressions given in Eqs. (2.1) and (2.2) guarantees that the estimators 
are unbiased (even in the presence of instrumental noise) if the Cp '° * (where Y, Y' are 



^We stress here that, strictly speaking, it is the adopted APS estimator to be optimal. In particular we 
have considered the QML method which is called optimal [35, 36] since it is unbiased and minimum variance 
(it provides the smallest error bars provided by the Fisher-Cramer- Rao inequality). Sec also Subsection 2.2 
and 2.3 and reference therein. 



any combination of T, E, B) are obtained through an unbiased angular power spectrum 
estimator ^ . 

The estimators Dj^ and Df^ have zero expectation value as can be readily seen by 
averaging Eqs. (2.1), (2.2) and using Eqs. (1.1), (1.2), (1.3), (1.4), (1.5). Therefore 

(dD = 0, (2.3) 

(^D = 0, (2.4) 

for every angle a. 

In order to find the most probable angle a compatible with the WMAP 7 year low 
resolution data, we adopt a Bayesian approach and maximize the log-likelihood defined as 

-21n£^(a) = ^^f '°'^M,f,^"'D,T'"'' , (2.5) 

w 

where X = TB or EB, D ' stands for Eq. (2.1) or (2.2) specialized to observed data (i.e. 
WMAP 7 year low resolution data in the following) and M^,^ = {Df D^). Note that in 
Eq. (2.5) we have already taken into account Eqs. (2.3), (2.4). 

2.2 Description of the simulations performed 

In this Section we describe the simulations we have performed in order to sample the 
likelihood defined in Eq. (2.5). 

In principle for each fixed angle a = a one has to simulate a large number of primordial 
maps, rotate them by a quantity 2(5 and add a random noise realization compatible with 
the considered data, as given in Eqs. (1.6), (1.7), (1.8). Then an angular power spectrum 
estimator has to be applied to each of the simulated maps to obtain the corresponding set of 
Cg '" *, where Y', y = T, E, B. Sue estimates are therefore used to build the matrix M^, 
for the angle a. Once M^, is inverted and contracted with JJ ' (of course computed for 
the same angle a), from Eq. (2.5) one obtains C{a). In order to sample £, this procedure 
has to be repeated for every angle a belonging to the range of interest [amin-, c^max]- 

In practice, since the likelihood given in Eq. (2.5) is by construction a periodic function 
of a with a period of 90°, one has to perform a MonteCarlo (MC) of Ngims simulations for 
each value of a between —45° and 45°. If for example the step (Aq) of such a sampling is 
1 degree, in principle one has to perform 91 MC of Ngims simulations each. This is not a 
short task since the computation time of one of those MC with Ngims = 10000 is around 4 
hours on 64 processors using an optimal estimator at low resolution (see [32, 33, 34, 35, 36] 
for further details on the implementation of the APS estimator we adopt). 

In fact, there is a shortcut to such a strategy that allows one to sample the like- 
lihood performing just one MC independently of the choice of the step Act. Feeding 



^This is not a trivial property. For example, note that the following non-linear estimator 

sm(2Q) = C^ ' / \{C^ ) +{C( • ) J 

is always biased, even for a noiseless experiment. 



Eqs. (1.6), (1.7), (1.8) into Eq. (2.1) and (2.2), after some algebra one obtains 

Df ^ = (Cf'"^ + <5C;^'"^) cos(2a) - (Cj'"^ + 5C;^'"^) sin(2a) , (2.6) 

j^EB ^ cf •"^cos2(2a) + <5C;^'"^ - ^ (cf '"^ + JC,"^'"^ + dCf'"'') sin(4a) , (2.7) 

where C^ '" is the angular power spectrum of the correlation between Y' and the noise 
of Y and 5C^ '" is the angular power spectrum of the residual noise leftover after the 

, . 1 . s;rinY',nY rinV ,nY ,nY',nY .,, ,nY',nY u ■ j-i, ■ i,- 

noise bias removal, i.e. oG^ = C^ — b^ with o^ being the noise bias 

estimate. Eqs. (2. 6), (2. 7) indicate clearly that the matrix M^,^ can be built starting from 
the random extractions of C^ '" and JC" '" which do not depend on a. 

2.3 Power spectra estimates from WMAP 7 data 

In this Section we describe the data set that we have considered. We use the temperature 
ILC map smoothed at 9.8 degrees and reconstructed at HealPix^ [37] resolution^ ^side = 16, 
the foreground cleaned low resolution maps and the noise covariance matrix in {Q, U) 
publicly available at the LAMBDA website ^ for the frequency channels Ka, Q and V 
as considered by Larson et al. in [38] for the low (. analysis of the WMAP data. These 
frequency channels have been co-added as follows [39] 

mtot = CtotiC^lrriKa + Cq^vtiq + Cy^my) , (2.8) 

where mj, Ci are the polarization maps and covariances (for i = Ka,Q and V) and 

Ctol = Qi + Cq' + Cy' . (2.9) 

This polarization data set has been extended to temperature considering the ILC map. We 
have added to the temperature map a random noise realization with variance of l/iK^ as 
suggested in [40]. Consistently, the noise covariance matrix for TT is taken to be diagonal 
with variance equal to IfiK^. 

We have performed Monte-Carlo simulations in order to properly sample the likelihood 
given in Eq. (2.5). A set of 10000 CMB plus noise sky realizations has been generated: 
the signal extracted from the WMAP 7 years best fit model, the noise through a Cholesky 
decomposition of the noise covariance matrix. We have then computed the APS for each 
simulation by means of our implementation of the QML estimator, namely BolPol [34], 
keeping record not only of the estimated Cj ^ but also of signal and noise only extractions. 
By saving the latter, one may obtain C^ '" and 5C^ '" given in Eqs. (2. 6), (2. 7). 

We considered two different Galactic masks: one for temperature and one for linear 
polarization, as suggested by the WMAP team [38]. Best-fit monopole and dipole have 
been subtracted from the observed ILC map through the HealPix routine remove-dipole 
[37]. 



*http://healpix.jpl. nasa.gov/ 

^For the reader who is not famihar with HEALPix convention, N^ide ~ 16 corresponds to maps of 3072 
pixels on the whole sky. 

® http://lambda.gsfc. nasa.gov/ 
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Figure 1: Left column of panels: likelihood for D^^. Middle column of panels: likelihood for 
D^^ . Right column of panels: joint likelihood for D^^ and D^^ . First row A£ = 2 — 6. Second 
row A£ — 2 — 23. Third row A£ — 2 — 47. Each panel represents the likelihood (normalized to unity 
at the peak) vs a, expressed in deg. See also the text. 



3. Results from WMAP 7 year 



3.1 Likelihood for a 



In Fig. 1 we show the likelihood distribution we obtain by sampling Eq. (2.5) for X=TB 
(left column), X=EB (middle column). Right column of Fig. 1 shows the joint likelihood 
of the two estimators computed as 



-21n/:^°^"*(a) = -2ln£^^{a) - 2 In/: 



-EB, 



a) 



(3.1) 



The best constraint we obtain for the birefringence angle comes from the joint likelihood 
for the largest range of multipoles (i.e. Ai = 2 — 47 corresponding to the lower-right panel 
of Fig. 1) and it is given by a = -L6° ± 1.7° (3.4°) at 68% (95%) C.L.. Considering 



Ai = 2 — 23 we obtain a 



^o+2.6° 
-2.5° 



i C.L. and a 



-3.0 



0-1-6.9° 
-4.7° 



at 



^0 C.L.. This 



is the same multipole range considered by the WMAP team in [17] where with a pixel 
based likelihood analysis they obtain a'^^^^'^y' = -3.8° ± 5.2° at 68% C.L.. 

We note that such results are compatible well within la and moreover we argue that 
the linear estimators of Eqs. (2.1) and (2.2) provide tighter constraints (at least when 
coupled with an optimal estimators for APS) with respect to the technique adopted in [17]. 

All errors quoted above and in what follows are purely statistical. We notice, inciden- 
tally, that the WMAP team quotes an instrumental systematic error on the polarization 
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10 20 30 40 

Figure 2: Spectrum of a vs £ obtained from Dj^ . Tlie dot represents tlie peak; position and the 
error bars provides the la C.L.. See also the text. 

angle at ±1.5° of the nominal orientation [17, 41, 42]. This further uncertainty contributes 
systematically to our final estimates. 

3.2 Angular spectrum for a 

As performed in [22] we specialize Eq. (2.5) to a given multipole range [^mim^max]) 

t-niax 

-2 In 4_,w) («) = E E Df-'^M§^-'Df^'^ , (3.2) 

where we have made explicit the dependence of irnm, ^max with only £' being summed over 
the entire available multipole range (i.e. from 2 to 47). This allows one to compute the 
spectrum of a, i.e. a vs i. In Fig. 2 we show the spectrum we obtain when we use Eq. (3.2) 
for X = TB. The error bars are obtained taking the la C.L.. 

Note that for X = EB at £ > 8 — 9 pure noise fluctuations overtake the signal-noise 
ones. Therefore, for such range of multipoles the sampled likelihoods might not be well 
featured. This is not the case for X = TB where this transition happens for ^ ^ 40 (see 
Appendix A for further details). This is the reason why we show the spectrum for TB 
in Fig. 2 but not the analogous for EB. The data for EB are taken into account only in 
combination with TB when building the joint likelihood, see Fig. 1. However, when the 
lowest multipoles are included in the analysis any possible noisy feature on the C^^ due 
to higher multipoles is subdominant, as shown by the middle column of Fig. 1. 

4. Planck forecasts 

In this Section we take into account the white noise level for 143 GHz channel of the 
Planck mission [43] launched into space on the 14th of May of 2009. As in [32, 33], we 
consider the nominal sensitivity of the Planck 143 GHz channel, taken as representative 
of the results which can be obtained after the foreground cleaning from various frequency 
channels. The 143 GHz channel has an angular resolution of 7.1' (FWHM) and an average 
sensitivity of Q ^K (11.4 ^i^) per pixel - a square whose side is the FWHM size of the 
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Figure 3: Joint likelihood of a for Ai = 2 — 16 in the left panel, Ai = 2 — 23 in the middle panel 
and Ai = 2 — 47 in the right panel. Red is for the WMAP data and blue for the Planck forecast. 
See also the text. 



beam - in temperature (polarization), after 2 full sky surveys. We assume uniform uncor- 
related instrumental noise and we build the corresponding diagonal covariance matrix for 
temperature and polarization, from which, through Cholesky decomposition we are able to 
extract noise realizations. For this low noise level we apply the same procedure adopted 
for the Monte-Carlo simulations in Subsection 2.3. 

In Fig. 3 we show the Planck forecast for the joint likelihood (plotted in blue). In 
the same panels we give also the WMAP joint likelihood for comparison (plotted in red). 
Left panel of Fig. 3 is for Ai = 2 — 6, middle panel for Ai = 2 — 23, and right panel for 
Ai = 2 — 47. The Planck likelihoods are computed as if Planck would have observed the 
same data as WMAP but with its own noise amplitude. We estimate an improvement by 
a factor 16.3 on the standard deviations. Of course the improvements shown in Fig. 3 are 
based only on the statistical noise and do not consider the systematic errors that are not 
publicly available now. However, following the analysis performed in [23] we can infer that 
the accuracy with which the angles and calibration of the Planck detectors must be known 
for the stated Planck forecasts to be dominant in the error budget, is Aa^^^^ < 0.1°. 



5. Conclusions 



We have performed an analysis of the WMAP 7 low resolution data in search for a birefrin- 
gence effect in the CMB. The large angular scale considered here probes the reionization 
era, as opposed to other higher resolution analyses that are more sensitive to last scattering 
epoch. The combination of an optimal estimator for the CMB angular power spectra and 
a linear estimator, Eq. (2.1) and (2.2), provides constraints in the range Ai = 2 — 23 
which are tighter by a factor ~ 2 than WMAP's own analysis, the only other results 
available in the literature at these large angular scales. We report our best constraint 
a = -1.6° ± 1.7° (3.4°) at 68% (95%) C.L. considering Ai = 2 - 47. We also provide, for 
the first time here, an angular spectrum of the birefringence angle a over the multipole 
range i < 50, detecting no significant deviation from the hypothesis of a scale independent 
a. We have, finally, forecasted the improvement expected from Planck data, finding that 
it should outnumber WMAP by a factor above 16 in terms of standard deviations. 
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A. Montecarlo simulations 

In this appendix we give some details about the angular dependence of the fluctuations 
entering the definition of our estimators, i.e. Eqs. (2.1) and (2.2). In Fig. 4 we present 
the averages and the standard deviation of the MC for the WMAP 7 year low resolution 
data. In particular, for each multipole and for each spectrum, we provide information 
about 5C" '" (left vertical bars), C^ '" (right vertical bars) and their sum (middle 
vertical bars). The averages are zero because the pure signal contribution does not enter 
the definition of the estimators and has been removed. It is clear that we have two kinds 
of behavior in Fig. 4: for X = TB (that means considering the panels for TB and TE) 
we have the signal-noise APS above the noise-noise APS up to £ ~ 40 (for the WMAP 
case) whereas for X = EB (that means considering the panels for EE, BB and EB) the 
signal-noise APS are of the same order as the noise-noise APS up to £ ~ 10 (except for BB 
whose signal is zero by construction). Beyond this scale the noise- noise APS is dominant. 
For the Planck case, see Fig. 5, the behavior is analogous to what already described 
for the WMAP case with the following differences: for X = EB the signal-to-noise ratio is 
improved for the lowest multipoles but the transition to noise dominated regime happens 
almost at the same critical multipole, i.e. ^ ~ 10 — 11 (since the primordial signal vanishes 
at that scale), and for X = TB the signal- noise APS are dominant for all the range of 
multipoles that are accessible at the considered resolution. Therefore, for the Planck data, 
when considering X = TB it might be worth to perform this analysis at higher resolution, 
e.g. Nside = 32. 
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